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was  based  on  physical  equations  for  a homogeneous  atmosphere  that  were 
modified  by  empirical  observations  made  under  a broad  range  of  meteorological 
conditions.  With  the  model  it  is  possible  to  estimate  the  visible  and  near- 
infrared  spectral  absorption  and  scattering  phenomena  that  result  from  atmo- 
spheric constituents,  from  surface  synoptic  weather  observations.  The  atmo- 
spheric model  was  coded  into  a computer  program  called  SCAT3. 

This  document  is  the  fourth  volume  documenting  the  radiometric  data  collec- 
tion program  that  culminated  in  the  radiometric  atmospheric  model,  SCAT3 . The 
objective  of  this  volume  is  to  acquaint  the  user  with  the  proper  input  sequences 
and  options  when  SCAT 3 is  run  on  an  IBM  370  computer.  Volumes,  1,  2,  and  3 
contain  the  details  of  data  collection,  analysis,  and  derivation  of  equations 
that  were  coded  into  the  program  to  reconstruct  atmospheric  radiance,  trans- 
mittance, and  ground  irradiance. 
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SECTION  1 

GENERAL  DESCRIPTION  OF  SCAT 3 PROGRAM 

1 . 1 INTRODUCTION 

SCAT3  is  a special  version  of  Program  SPECTR,  a general  purpose  program 
for  manipulating  optical,  photographic,  and  electro-optical  elements 
on  a wavelength-by-wavelength  basis.  Each  element  is  represented  by  an 
array  of  values  representing,  for  instance,  transmittance,  reflectance, 
or  relative  sensitivity  for  a 5 nm  portion  of  the  spectrum  (300-1200  nm) . 
These  elements  may  be  added,  multiplied,  integrated,  etc,  to  model  the 
response  of  a system.  A disk-based  library  of  elements  provides  flexi- 
bility in  modeling  a wide  variety  of  components  and  systems.  In  addition 
to  handling  array  values,  the  program  may  also  use  temporary  constants. 

SCAT3  contains  an  atmospheric  simulation  model  based  on  the  collection 
program  and  data  analysis  discussed  in  detail  in  Volumes  1,  2,  and  3. 

The  atmospheric  model  is  contained  in  a series  of  subroutines,  commanded 
by  an  operator  command  function,  AMOS.  The  output  of  the  AMOS  function 
is  a series  of  spectral  elements  describing  atmospheric  path  transmittance 
and  radiance,  various  forms  of  ground-level  irradiances,  and  several  other 
spectral  radiometric  quantities  that  can  be  used  in  subsequent  computations. 
User  function  SYSINT  VFUN  is  available  for  performing  a series  of  calcula- 
tions and  integrations  of  AMOS  elements, and  elements  that  define  a system- 
sensor response. 

1 . 2 OPERATING  COMMANDS 

The  SCAT3  Program,  in  its  basic  form,  supports  four  commands  which  in  turn 

support  the  functions  shown  in  Figure  1-1.  These  commands  are  summarized 

below  and  described  in  detail  in  the  following  four  sections  of  this  volume: 

a.  AMOS  generates  temporary  elements  for  atmospheric  computations, 
based  on  surface  weather,  solar  altitude,  etc. 
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b.  GENERATE  generates  temporary  elements  that  are  the  result  of 
computations;  e.g.,  special  functions  such  as  AMOS,  or  the 
user  function  SYSINT  VFUN.  The  command  GENERATE  SCAT3  AMOS 
is  equivalent  to  AMOS. 

c.  I/O  (Input/Output) is  used  to  read  in  or  print  out  spectral 
elements . 

d.  FILE  is  used  for  disk  file  maintenance.  It  is  normally  used 
to  save  spectral  elements  in  the  library  or  delete  elements 
when  they  are  no  longer  needed. 


Although  SCAT3  has  a multiplicity  of  optional  functions,  only  the  follow- 
ing four  basic  command  sequences  are  used: 

a.  GENERATE  SC AT 3 AMOS, or  simply  AMOS 

b.  GENERATE  SYSINT  VFUN 

c.  I/O  READIN 

d.  FILE  SAVE  or  FILE  DELETE 


The  various  other  options  are  not  used  in  routine  operation  of  the  SCAT3 
program  and  have  therefore  been  grouped  for  further  discussion  in  Section  6. 


1.3  EXAMPLE  OF  A TYPICAL  PROGRAM  RUN 


To  set  up  program  SCAT3,  a series  of  command  and  associated  data  cards  aie 
used.  Each  command  card  (e.g.,  GENERATE  ...)  must  begin  in  column  1 of 
the  card.  If  the  particular  function  does  not  require  a data  card,  subse- 
quent command  cards  for  the  same  function  may  have  a blank  in  column  1 
instead  of  the  command  name.  An  80-column  layout  sheet  for  the  command  and 
data  cards  for  a typical  run  are  shown  in  Figure  1-2.  The  purpose  of  each 
command  sequence  in  the  example  in  Figure  1-2  is  summarized  below.  This  se- 
quence exercises  all  the  required  commands  in  the  normal  operation  of  SCAT3. 

(1)  Read-in  a new  detector  spectral  sensitivity  element  to  be 
stored  by  the  name  DET78. 

(2)  Peak-normalize  DET78,  and  store  the  normalized  spectral 
element,  by  the  name  DET78,  in  the  disk  file. 
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(3)  Generate  the  atmospheric  and  ground  level  irradiance 
spectral  elements  corresponding  to  the  following: 
solar  altitude  of  40  degrees,  surface  temperature 

of  15  degrees  Celsius,  relative  humidity  of  60%, 
forest  surround,  and  all  other  parameters  having 
default*  values  for  the  date  4/1/78. 

(4)  Spectral,  integrate  the  elements  DH,  T,  and  SK,  etc, 
with  the  default  spectral  elements  for  detector,  lens, 
and  scene  reflectance. 

(5)  Repeat  the  spectral  integrations,  substituting  the 
normalized  detector  array  PKDET78  for  its  default 
element.  Multiply  the  results  by  the  factor  1.375. 

The  lens  and  scene  reflectance  elements  remain  un- 
changed. 

(6)  Delete  the  element  DET77  from  the  disk  library. 


If  no  input  specification  is  given,  the  program  will  assume  certain 
values  of  parameters  and  spectral  element  arrays. 
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DESCRIPTION  OF  THE  COMMAND  AMOS 

2.1  GENERAL  DESCRIPTION 

The  AMOS  command  allows  the  user  to  generate  spectral  elements  that 
describe  the  atmosphere  based  on  estimated  surface  meteorology,  solar 
altitude,  acquisition  geometry,  and  other  factors.  These  spectral  elements 
can  be  used  in  subsequent  functions,  such  as  SYSINT  VFUN,  I/O  FILE,  etc. 

2.2  INPUTS 

Each  run  of  AMOS  requires  a command  card  that  defines  print  options,  and  a 
data  card.  The  data  card  provides  the  values  for  the  input  variables  that 
are  needed  to  generate  the  spectral  atmospheric  and  ground  irradiance  elements. 
All  options  and  values  have  defaults  except  for  solar  altitude  and  weather/ 
transmittance  estimates.  The  input  values  are  checked  for  reasonableness, 
and  warnings  or  fatal  error-condition  messages  may  be  written,  depending  on 
the  seriousness  of  the  error. 

The  user  may  specify  the  amount  of  printout  from  an  AMOS  run.  The  choices 
are:  FULL,  which  prints  all  8 spectral-array  elements;  PARTIAL,  which 
prints  only  atmospheric  path  transmittance  (T) , atmospheric  path  radiance 
(SK) , and  horizontal  daylight  irradiance  (DH) ; and  NONE,  which  prints  no 
array  elements.  The  input  and  default  variables  are  printed  regardless  of 
the  option  selected. 

2.2.1  Command  Card 

The  three  valid  AMOS  commands  that  serve  to  illustrate  the  print  options 
(see  Figure  2-1)  are: 

a.  GENERATE  SCAT3  AMOS  NONE,  or  AMOS  NONE 

b.  GENERATE  SCAT 3 AMOS  PARTIAL,  or  AMOS  PARTIAL 

c.  GENERATE  SCAT 3 AMOS  FULL,  or  AMOS  FULL. 
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2.2.2  Data  Card 


Each  AMOS  command  card  must  be  followed  by  a single  data  card  having  the 
format  and  definition  of  input  variables  shown  in  Tables  2-1  and  2-2.  When 
a run  is  set  up,  the  user  must  be  sure  that  any  values  that  he  lets  default 
match  the  desired  run  conditions.  The  program  can  be  operated  using  one  of 
three  types  of  input  variables. 

Surface  Weather  for  which  the  inputs  are  surface  relative  humidity, 
temperature  in  degrees  centigrade,  and  pressure  corrected  to  sea 
level , in  mi  1 1 ibars . 

T550  for  which  the  single  input  is  an  estimate  of  the  vertical-path 
atmospheric  transmittance  at  550  nm  wavelength,  and 

Precipitable  Moisture  for  which  the  simple  input  is  the  integrated 
partial  pressure  of  water  in  a vertical  column  of  atmosphere,  ex- 
pressed in  centimeters. 

These  input  options  are  mutually  exclusive,  and  if  input  variables  are 
specified  for  more  than  one  mode,  error  messages  will  be  generated  and  the 
AMOS  run  will  possibly  be  terminated.  Valid  data  cards  that  illustrate  use 
of  the  AMOS  atmospheric  options  are  shown  in  Figure  2-2  . 
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FORMATS  AND  DEFINITIONS  OF  DATA  CARD  PARAMETERS 
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The  decimal  reflectance  at  650  mm,  which  is  multiplied  by  a spectral  array  selected  bv 
ITYP  that  is  normalized  at  650  mm. 


2.3  OUTPUT 


The  output  of  SCAT3  is  printed  copy  only.  Each  execution  of  the  program, 
which  results  from  a correctly  input  command  and  a data  card,  generates  the 
output  described  below. 


2.3.1  Data  Card  Variables 

The  listing  of  variables  is  divided  into  six  groups  accordings  to  the  way 
they  operate  in  the  program  and  as  shown  in  Figure  2-3.  The  six  variable 


groups  are: 

1.  Trigonometric  Data  (in  degrees) 

SA: 

Solar  altitude 

CATS: 

Observer-sun  angle 

ZA: 

Observation  xenith  angle 

LAT: 

Station  latitude 
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Figure  2-3.  Example  of  SCAT3  Output  and  Listing  of  Data-Card- 
Specified  Variables 


2.  Surface  Weather  Data 


TEMP: 

RH: 

SL  PRES: 
DUST: 


Temperature  (*C) 

Relative  humidity 

Station  pressure  corrected  to  sea- level  (mb) 
Surface  dust  concentration  (part/cm3) 


3 . Alternative  Scaling  Parameters 

T550:  Vertical  atmospheric  transmittance  at  550  nm 
02:  Ozone  scale  height  (cm) 

PM:  Total  vertical  column  of  precipitable  moisture; 

sea- level  to  space  (cm) 

4 . Elevation  Data  (km) 

ELEV:  Station  elevation  above  sea- level 

VA:  Observer's  altitude  above  station  elevation 

5.  Surround  Reflectance  Parameters 

ALB  TYPE:  Description  of  station  surrounding 

(i.e.,  vegetation,  sand,  snow,  industrial) 

ALB  REF:  Surround  reflectance  at  650  nm 

6 . Atmosphe r ic  Tra nsmittance  P arameters 

KR:  Rayleigh  scattering  coefficient 
KA:  Aerosol  scattering  coefficient 
C:  Exponent  for  aerosol  scattering  term 

W:  Precipitable  moisture  between  station  and 

observer’s  elevation  (cm) 

2.3.2  Spectral  Arrays 


In  addition  to  the  data  card  variables,  the  spectral  arrays  and  their  numeric 
codes  that  were  generated  as  a result  of  these  input  data  card  variables,  are 
also  printed. 

These  arrays  can  be  recalled  from  temporary  program  data  files  to  operate  on 
the  spectral  data  in  subsequent  calculations.  The  arrays  are  not  permanently 
stored  and  are  written  over  upon  the  execution  of  a subsequent  GENERATE  command. 
The  eight  spectral  arrays  grouped  by  type  are  listed  below. 
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Atmospheric  Path  Transmittance  (350-11S0  nm) 


T:  Oblique  spectral  transmittance  (scattering  plus  absorption) 

from  the  observer’s  altitude  to  ground,  adjusted  for  precipitable 
moisture,  dust,  path  length,  and  terrain  elevation. 

TV:  Vertical  spectral  transmittance  (scattering  plus  absorption) 

from  the  observer's  altitude  to  ground,  adjusted  for  precipitable 
moisture,  dust,  and  terrain  elevation. 

Atmospheric  Path  Radiance  (350-1 ISO  nm) 

SK:  Spectral  atmospheric  path  radiance,  at  the  specified  angles, 
adjusted  for  observer's  altitude,  terrain  elevation,  surround 
albedo,  precipitable  moisture, (watts/m2/steradian/5  nm) . 

Terrain  Reflectance  (350-1150  nm) 

ALB:  Spectral  reflectance  of  the  station  surround  based  on  the  albedo 

type  and  normalizing  reflectance  that  is  input. 

Ground  Level  Irradiance  (350-1150  nm) 

DH:  Daylight  irradiance  on  a horizontal  surface  adjusted  for  solar 
altitude,  atmospheric  transmittance,  and  surround  albedo 
(watt/m2/5  nm) . 

SH:  Skylight  irradiance  on  a horizontal  surface  from  the  full  hemisphere 
minus  direct  sunlight  adjusted  for  solar  altitude,  atmospheric 
transmittance,  and  surround  albedo  (watts/m2/S  nm) . 

DV:  Daylight  irradiance  on  a vertical  surface  perpendicular  to  the 

solar  azimuth  adjusted  for  solar  altitude,  atmospheric  transmittance 
and  surround  albedo  (watts/m2/5  nm) . 

SVB:  Skylight  irradiance  on  a vertical  surface  perpendicular  to  and 

180®  away  from  the  sun  azimuth  adjusted  for  solar  altitude,  atmo- 
spheric transmittance,  and  surround  albedo  (watts/m2/5  nm) . 

Figure  2-4  shows  a typical  spectral  array  generated  by  the  AMOS  subroutine. 
Although  only  three  such  arrays  are  printed  under  the  PARTIAL  print  option, 
all  eight  arrays  are  internally  generated  and  available  for  further  computer 
operations. 
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2.3.3  Error  Messages 


If  a parameter  specified  on  the  Data  Card  is  out  of  range  or  a required 
option  is  omitted,  an  appropriate  error  message(s)  will  precede  the  reprint 
of  the  command  card  values.  Figure  2-5  shows  an  output  page  that  lists 
several  possible  error  messages.  This  is  not  a complete  list  of  all  possible 
messages  but  it  illustrates  their  format.  A complete  listing  of  possible 
error  messages  generated  by  subroutine  AMOS  is  contained  in  Appendix  B. 
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NOTE:  The  improperly  specified  data  card  below  was  responsible 
for  the  error  messages  above. 


SECTION  3 

DESCRIPTION  OF  COMMAND:  GENERATE  SYSINT  VFUN 

3.1  COMMAND  FORMAT 

A special  user  function  SYSINT  VFUN  has  been  provided  that  calculates 
spectrally  integrated  ground- level  irradiances,  atmospheric  path  radiance, 
and  effective  atmospheric  transmittance  for  a generalized  system.  This 
function  is  used  in  conjunction  with  an  AMOS  run,  which  must  precede  it, 
however.  By  overriding  the  default  elements  and  the  system  constant,  the 
computations  can  be  performed  for  a wide  variety  of  system  components. 

The  default  elements  are  listed  in  Appendix  C.  Figure  1-2,  items  3,  4,  and 
5 illustrate  the  input  for  a typical  run  of  this  function  along  with  its 
associated  AMOS  run. 

3.2  INPUTS 

A typical  command  is:  GENERATE  SYSINT  VFUN  XXX  YYY  ZZZ . The  variables  XXX, 
YYY,  and  ZZZ  are  element  names  for  the  detector,  optics,  and  ground  reflec- 
tance, respectively.  These  elements  can  be  extracted  from  the  disk  library 
or  entered  on  cards  via  a preceding  I/O  READIN  step.  If  any  of  these  elements 
are  defined  short  of  the  range  400  to  1100  nm,  the  integration  limits  are 
adjusted  accordingly.  The  element  defining  the  detector  should  be  peak- 
normalized  to  the  value  of  1.  The  user  can  set  all  three  elements  to  their 
default  values  by  replacing  XXX  with  STD  and  omitting  YYY  and  ZZZ.  Any  single 
element  can  be  set  to  its  default  value  by  specifying  in  DEF  in  place  of  the 
element  name.  This  is  shown  in  Figure  3-1.  The  variables,  which  are  specified 
as  XXX,  YYY,  or  ZZZ  use  the  input  codes  defined  in  the  table  below. 

An  optional  units-conversion  multiplier  may  be  specified  as  K = F.FF.  The 
printed  values  of  output  ERS  and  EH  will  be  multiplied  by  F.FF.  If  no  value 
is  specified,  1.00  is  assumed. 


■ 


Element 

Names 


XXX  (for 
Detector) 


Codes 


STD 


DEF 


Any  valid 
element  name 
(i.e,835N) 


I i 


YYY  (for 
Optics) : 


ZZZ  (for 

Ground 

Reflector) 


DEF 


Any  valid 
element  name 
(i . e. ,424L) 


DEF 


Any  valid 
element  name 
(i.e. ,179RF) 


! 


Definition 


The  code  that  permits  all  three  spectral 
elements  to  assume  their  default  values. 

(XXX  becomes  899N,  YYY  becomes  457L,  and 
ZZZ  becomes  825F.) 

The  code  that  only  permits  XXX  to  default  to 
899N. 


A spectral  array  defined  in  the  library  of 
elements  that  represents  the  system  spectral 
response.  The  values  should  be  peak-normalized. 


The  code  that  only  permits  XXX  to  default  to 
457L. 

A spectral  array  defined  in  the  library  of 
elements  representing  the  system  spectral 
optical  transmittance. 


The  code  that  only  permits  ZZZ  to  default  to 
825F. 


A spectral  array  defined  in  the  library  of 
elements  representing  the  spectral  reflectance 
of  a ground  object. 
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3.3  OUTPUTS 


Each  use  of  the  SYSINT  VFUN  Command  generates  one  page  of  output.  An 
example  of  this  output  is  shown  in  Figure  3-2.  The  format  for  this  output 
is  summarized  below. 


Line(s) 

1 

2-5 


6 

7 


Data  Provided 

Array  element  names  for  detector  (SD) , optics  (OPT), 
and  reflectance  on  the  ground  (REF) . 

The  integrated  effective  value  for  the  rudiance  (ERS) 
assumed  reflector  (REF)  under  the  various  forms  of 
ground  level  irradiance  and  effective  path  transmittance 
(TEFF)  relative  to  each  irradiance. 

The  integrated  effective  atmospheric  path  radiance  x n (l.ll) 
The  assumed  or  read-in  value  of  K. 


If  one  or  more  of  the  spectral  elements  is  not  defined  over  the  range  400  to 
1100  nm,  the  integration  limits  are  changed.  If  this  occurs,  a warning  message 
e.g.,  UPPER  INTEGRATION  LIMIT  LOWERED  TO  1000  nm,  is  printed. 

The  two  following  error  indications  may  possibly  be  encountered. 

(1)  NO  AMOS  ELEMENTS  LOCATED  . . . which  indicates  that  the  program 
was  unable  to  locate  elements  T,  SK,  and  DH,  etc,  from  a prior 
AMOS  run  or  from  the  disk  library. 

(2)  UNABLE  TO  LOCATE  XXX,  which  indicates  that  the  spectral  array 
element  XXX  was  not  found  in  the  disk  library  or  from  a pre- 
vious step. 

3.4  SUMMARY  OF  COMPUTATIONS 


The  following  equations  are  used  in  SYSINT  VFUN 
a.  Effective  Integrated  Irradiances: 

[1100 

J H(X)  • Rs(X)  • T(X) 
400 


k(X) 


SD(X)  dX 


] 
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Figure  3-2.  Example  of  SCATS  Output:  Integrated  Values  Generated 
by  SYS 1 NT  VFUN  Command 


b.  Atmospheric  Path  Radiance: 


EH 


it  • K 


1100 

/ "hW 


400 


L(A) 


Sd(A)  dA 


c.  Effective  Atmospheric  Path  Transmittance: 


TEFF 


1100 

f H(A) 
400 


1100 


/ 


H(A) 


400 


RS(A)  • T(A)  • L(A)  • Sd(A)  dA 


RS(A)  • L ( A ) • Sd(A)  dA 


In  the  above  equations,  the  integrand  variables  are  defined  as  follows: 


T(A)  is  the  spectral  atmospheric  path  transmittance 
(AMOS  array  name  is  T). 

N^(A)  is  the  spectral  atmospheric  path  radiance  in  watts/m2/S  nm 
steradian  (AMOS  array  name  is  SK). 

H(A)  is  either  skylight  or  daylight  irradiance  on  either  a horizontal 
or  a vertical  surface  in  watts/m2/5  nm  (AMOS  array  names  are  DH, 
DV,  SH,  SVB) . 


In  addition,  the  integrand  variables  from  the  disk  library  or  a preceding  I/O 
READIN  step  are  defined  as  follows: 

Sq(A)  is  the  spectral  sensitivity  of  the  Film/Detector  (peak-normalized) 

L(A)  is  the  spectral  transmittance  of  the  optics,  and 

Rg(A)  is  the  spectral  ground  reflectance. 
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SECTION  4 


DESCRIPTION  OF  COMMAND  I/O  READ IN 


4.1  GENERAL  DESCRIPTION 

This  command  is  used  to  read  any  element  array  into  the  program  for  use  in 
computations  or  for  permanent  storage  in  the  disk  library.  The  general  form 
of  the  command  is: 

I/O  READIN  (option) 

The  possible  options  are: 

NOLIST  which  allows  the  spectral  array  to  be  read  into 
a temporary  file,  but  returns  no  printed  listing 
of  the  array. 

SAVE  which  allows  the  spectral  array  to  be  read  in  and 
stored  as  a permanent  array  by  the  element  name 
given  on  the  header  card  (equivalent  to  a FILE  SAVE 
Command,  see  section  5),  or, 

"Blank"  which  allows  the  spectral  array  to  be  read  into  a 
temporary  file  and  results  in  a printed  listing  of 
the  array. 

Any  element  of  the  same  name  that  is  already  in  the  library  is  printed  out  and 
replaced.  The  response  SAVED  ON  DISK  indicates  satisfactory  completion.  Follow- 
ing the  command  card  are  a header  card  and  data  cards  that  contain  the  spectral 
values . 

4.2  INPUTS 

Once  the  I/O  READIN  command  has  been  specified,  a number  of  spectral  arrays 
can  be  read  in.  Each  array,  or  element,  is  read  in  using  a header  card  and  a 
series  of  data  cards  in  the  formats  described  below. 
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HEADER  CARDS 


Column 

1-8 

10-12 

20-22 

29-32 

40-49 

73-80 

Column 

1-70 

73-76 

79-80 


Format 


Description 


A8  Name  of  element  in  4 to  8 alphanumeric  characters, 

at  least  one  of  which  must  be  alphabetic. 


A3  One  of  the  following  four  transformations  must  be 

used  to  indicate  condition  of  the  data. 


VAL: 

In 

correct  form 

LOG: 

In 

logarithmic  form 

LIN: 

In  log-inverse  form 

INV: 

In 

inverse  form 

13  The  wavelength  of  first  data  point. 


14  (right  The  wavelength  of  last  data  point, 
justified) 

F10.4  The  constant  that  is  multiplied  times  every  data 
point  after  conversion  to  the  correct  form.  If 
left  blank,  default  is  1.0. 


A8  (right  The  date  when  the  deck  was  made.  It  must  always 
justified)  be  present  and  never  use  less  than  columns  75-80 
(ex:  1/1/78). 


DATA  CARDS 


Format 


Description 


10E7.2  Lists  spectral  data  starting  with  the  value  for 

the  first  wavelength  specified  on  the  header  card 
and  continuing  on  other  cards  until  the  value  for 
the  last  wavelength  has  been  reached. 

A4  Identification. 

12  (right  The  order  number  of  the  card, 
justified) 
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SECTION  5 

DESCRIPTION  OF  THE  COMMAND  FILE 


5.1  GENERAL  DESCRIPTION 

This  command  is  used  to  maintain  the  disk  library  of  elements.  Each  card 
must  contain  one  complete  command  to  store  or  delete  a spectral  array  on  the 
disk  file. 

The  first  card  must  contain  the  word  FILE  beginning  in  column  1.  The  cards 
immediately  following  the  first  card  need  not  contain  the  word  FILE,  but  in- 
stead can  be  blank  in  column  1.  Each  command  item  (SAVE,  DELETE,  NAME,  etc) 
must  be  separated  from  the  next  item  at  least  one  blank.  Blanks  must  not  be 
used  within  a command  word. 

The  name  by  which  an  array  is  to  be  saved  in  a disk  file  may  contain  (1)  as 
many  as  eight  characters  and  (2)  any  combination  of  letters  and  numerals. 
However,  all  numerals  cannot  be  used. 

5.2  INPUTS 

Three  possible  FILE  commands  can  be  given,  each  of  which  must  be  struc- 
tured as  indicated  in  paragraph  5.1  above. 

5.2.1  FILE  SAVE  XXX  BY  YYY 

(In  the  following  paragraphs,  XXX  refers  to  the  name  of  an  array 
that  is  either  to  be  stored  by  the  new  name  YYY  or  deleted;  XXX 
and  YYY  can  be  identical.) 

This  command  performs  the  basic  task  of  placing  XXX  in  the  disk  file  and 
giving  it  the  element  name  YYY.  If  another  array  is  already  on  the  disk  by 
the  second  name,  the  save  operation  will  not  be  executed.  Instead,  both 
arrays  will  be  printed. 


5.2.2  FILE  SAVE  XXX  BY  YYY  U 


This  command  differs  only  from  that  of  5.2.1  in  that  it  replaces  existing 
elements  in  the  disk  file  (e.g.,  with  an  improved  estimate  of  the  same  lens 
transmittance).  After  the  array  XXX  is  stored  by  the  name  YYY,  both  the 
old  element  YYY  and  the  new  element  YYY  will  be  printed. 

5.2.3  FILE  DELETE  XXX 

The  purpose  of  this  command  is  to  erase  the  specific  element  XXX  from  the 
disk  file.  After  the  array  is  deleted,  the  element  XXX  is  printed  to  indi- 
cate what  has  been  erased. 
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SECTION  6 
OTHER  FUNCTIONS 

6.1  GENERAL  DESCRIPTION 

In  addition  to  the  commands  already  discussed,  SCAT3  also  supports  several 
other  command  functions  that  are  occasionally  used.  These  commands  are  given 
on  a single  card,  and  they  usually  generate  one  page  of  output.  The  element 
that  is  created  by  these  commands  is  often  used  in  subsequent  steps.  Two  of 
these  functions  are  discussed  below. 

6.2  I/O  FUNCTIONS 

In  addition  to  the  I/O  READIN  function  discussed  in  Section  4 of  this  report, 
function  I/O  PRINT  and  I/O  DUMP  are  also  available.  The  first  of  these,  I/O 
PRINT,  is  used  to  print  an  element  array  or  constant.  For  this  item,  the  com- 
mand is  I/O  PRINT  XXX  YYY  2Z2  . . . where  XXX,  YYY,  etc,  are  the  names  of  elements 
in  the  program,  or  in  the  disk  library.  The  user  may  specify  as  many  element 
names  as  will  fit  on  the  command  card,  but  each  name  must  be  separated  by  blanks. 
An  example  of  the  output  format  for  any  array  element  is  shown  in  Appendix  C« 

The  command  I/O  DUMP  is  used  by  the  programmer  who  maintains  the  disk  library; 
it  can  also  be  used  in  debugging  problems  within  the  program.  The  command 
I/O  DUMP  will  generate  listings  of  the  names  of  all  elements  in  the  disk  library, 
and  those  in  the  program's  core  data  table  (CDT) . 

6.3  GENERATE  FUNCTION 

In  addition  to  the  command  GENERATE  SYSINT  VFUN  described  in  Section  3 of  this 
report,  there  are  several  other  GENERATE  functions.  The  resulting  output  is 
printed  in  the  format  shown  in  Appendix  C for  an  array  element,  or  one  value 
per  page  for  a constant.  The  inputs  for  these  functions  are  of  two  types: 
single  operator  functions  and  double  operand  computations. 


6.3.1  Single  Operator  Functions 

These  commands  are  of  the  form:  GENERATE  XXX  ■ (function)  YYY,  where  the 
function  can  be  any  of  those  described  below. 

(1)  SUM.  Integrates  the  values  of  an  array  element  over  its 
spectral  range.  The  output  is  a constant. 

Example:  GENERATE  X ■ SUM  XYZ  where  XYZ  is  an  array  element 
and  X will  be  the  integral  constant. 

(2)  LOG.  Takes  the  log  (base  10)  of  a constant  or  of  each  value 
of  an  array. 

Example:  GENERATE  LOGX  « LOG  X. 

(3)  AREA  NORMALIZE.  Creates  a scaled  array  such  that  the  area 
under  the  curve  is  equal  to  1.0,  and  the  shape  is  the  same 
as  the  input  array. 

Example:  GENERATE  ANX  ■ AREANORMALIZE  X,  where  X is  a spec- 
tral array. 

(4)  PEAK  NORMALIZE.  Creates  a scaled  array  such  that  the  maximum 
value  is  1.0,  and  the  shape  is  the  same  as  the  input  array. 

Example:  GENERATE  PKK  « PEAKNORMALIZE  X,  where  X is  a spec- 
tral array. 

(5)  BBODY.  Generates  an  array  element  from  300-1200  nm  for  a black 
body  radiator  at  the  given  color  temperature  (°K) . 

Example:  GENERATE  BB  BBODY  3200. 

Note:  No  equals  sign  (-)  is  used. 

(6)  BANDPASS.  Generates  an  array  element  from  300-1200  nm  for  a band- 
pass filter.  At  wavelengths  from  LON  to  HIGH,  the  element  has  values 
of  1.  Elsewhere  it  has  values  of  xero. 

Example:  GENERATE  FILT  BANDPASS  (450)  (1025)  where  low  is  450  nm 
and  HIGH  is  1025  nm. 
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(7)  User  Function.  Program  SPECTR  provides  for  three  user-written 
functions  to  give  access  to  program  elements  and  utility  routines. 

The  user  function  is  a FORTRAN  program  with  the  name  SUBROUTINE  UFUNCT 
( , , ) . As  with  other  functions  in  the  program  SPECTR,  the  user  function 
must  be  called  by  an  appropriate  command  of  the  form: 

GENERATE  XXX*  CFUN  (option,  option,  . . .) 

SCAT  3,  for  example,  uses  a user  function**  routinely  to  perform  spectral 
integrations  of  arrays  generated  by  AMOS.  Within  the  constraints  of  subroutine 
name  and  command,  complete  flexibility  is  contained  in  this  user  option. 

Appendix  F.  reviews,  in  more  detail,  the  generation  of  a user  function  routine 
which  can  be  used  as  a guide  in  development  of  new  ones.  Reference  to  the  SCAT3 
program  listings  must  be  made  before  one  attempts  to  write  a user  function. 

6.3.2  Double  Operand  Functions 

The  GENERATE  Command  also  provides  for  double  operand  computations  such  as 
addition,  subtraction,  multiplication,  division,  and  exponentiation.  The  user 
may  perform  the  above  computations  on  constant  elements,  array  elements,  or 
numeric  values.  However, an  array  cannot  be  raised  to  an  array  power  or  divided 
by  zero. 

Examples  of  the  above  computations  are: 

GENERATE  A * B ♦ C 
X • Y * A 

Y * M * X 

Y * Y ♦ B 

PI  - 3.14159  x 1 
ANTILOG  « 10.  **  A 

The  word  GENERATE  may  be  omitted  after  the  first  GENERATE 
command  if  the  subsequent  commands  do  not  begin  in  column  1. 


A user-provided  variable  name  that  is  subject  to  the  restrictions 
discussed  in  Section  4.2.1,  Element  Name. 

See  Section  3 of  this  report. 
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APPENDIX  A 


SCAT3  TEST  CASES 


A series  of  IS  test  cases  was  constructed  and  run  to  utilize  all 
possible  input  options  and  variables  and  to  illustrate  the  output 
of  a correctly  assembled  input  deck.  Figure  A-l  is  an  annotated  copy 
of  the  input  commands  that  were  used,  and  Table  A-l  lists  the  re- 
spective 15  output  results. 
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TAV6S0.  131077100.  -.5  105.3 z.  Errors 


INTEGRAL  RESULTS  FOR  SCAT 3 TEST  CASES 

(K=l .0) 


10,000  ft.  altitude. 

Snow  surround  case  with  K • 0.983 


APPENDIX  B 

ERROR  MESSAGES  THAT  MIGHT  BE  GENERATED  BY  SCAT 3 


The  meaning  of  each  of  the  following  error  messages  is  explained  on  the 
line  below  the  message. 

1.  ALB  RF  XX. XX  OUT  RANGE  (0-1) 

Albedo  surround  reflectance  value  XX. XX  not  between  0 and  1. 

2.  ALB  TYPE  N INVALID 

Albedo  surround  type  N not  valid  (0-4). 

3.  YOU  MUST  SPECIFY  A TYP  AS  WELL  AS  AN  ALB  REF 

If  you  override  the  default  albedo  surround  reflectance  you  must  specify 
TYP  1-4. 

4.  AN  ELEVATION  OF  XX. X KM  SEEMS  TOO  HIGH 
Warning. 

5.  YOU  ARE  UNDERGROUND  XX. XX  KM 
A negative  observer  altitude. 

6.  LAT  XX. XX  INVALID 

Latitude  must  be  -90  to  +90.  For  southern  hemisphere  (negative  latitude) 
an  ozone  value  of  3.0mm  is  assumed. 

7.  SA  XX. X DEG  OUT  OF  RANGE  (-5  -90) 

The  solar  altitude,  after  any  low  solar  altitude  refraction  correction, 
exceeds  the  range  of  the  model. 

8.  RH  XX. XX  OUT  OF  RANGE  (0  - 1.00) 

Relative  humidity  out  of  range. 

9.  ZEN  ANGLE  XX. XX  OUT  OF  RANGE  (0  - 90) 

Observation  angle  not  between  0 and  90  degrees. 

10.  CATS  ANGLE  XX. X OUT  OF  RANGE  (0  - 180) 

11.  XX. X DEG  CELSIUS  OUT  OF  RANGE  (-/+50) 

Temperature  out  of  range. 

12.  SPECIFYING  PM  X.XXX  OVERRIDES  RH  XXX 
Warning. 

13.  YOU  MUST  SPECIFY  EITHER  RH,  PM,  or  T550. 
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14.  YOU  SPECIFIED  BOTH  A RAYLEIGH  ATMOSPHERE  & A PM  OF  X.XXX 

15.  YOU  SPECIFIED  BOTH  A RAYLEIGH  ATMOSPHERE  & A RH  OF  X.XXX 

16.  T550  NNNN  OVERRIDES  PM  XXX  &/OR  RH  X.XX 
Warning. 

17.  T550  X.XXX  OUT  OF  RANGE  (.2  - .815) 

18.  READ/FORMAT  ERROR  IN  AMOS  INPUT  CARD 
Error  return  from  FORTRAN  READ  (.  . . 

19.  ANGLE  IN  BEMAS  X.XX  GT  90  DEGREES 

Warning,  which  may  occur  with  negative  solar  altitude. 

20.  CONVERTED  X.XX  FT  TO  Y.YY  KM 
Advisory. 

21.  NO  INPUT  DATE  --  ASSUMED  TODAY 
Advisory. 


22.  INPUT  DATE  WRONG  - ASSUMED  YOU  MEANT 
Warning  when  day  and  month  interchanged. 

23.  DATE  NOT  USABLE 

24.  CREARH  -VE  PM  ...  PM,  KA,  K03,  TS  

Combination  of  T550,  ozone,  pressure,  temperature,  dust  etc,  give  a 
negative  value  for  calculated  precipitable  moisture.  This  error  would 
be  followed  by  a system  "log  of  negative  number"  error. 

25.  RH  X.XX  5 TEMP  YY.Y  DEG  C NOT  CONSISTENT  WITH  T550  Z.Z  IN  CREARH 
Calculated  a relative  humidity  over  100% . 

26.  DEFAULT  TEMPERATURE  INCREASED  TO  XX. X DEG  C 
Advisory. 

27.  DO  YOU  REALLY  WANT  OZONE  AT  THE  EQUATOR? 

Advisory  when  no  OZ  specified  and  LAT  is  zero. 

28.  PRESSURE  XX. XX  IN  HG  CONVERTED  TO  YY.Y  MB 

Advisory  when  units  converted  from  inches  of  mercury  to  millibars. 

29.  PRESSURE  X.XX  MB  OUT  OF  RANGE  (500  - 1200) 
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30.  AMOS  ELEMENT (S)  NOT  LOCATED 

Run  of  SYSINT  VFUN  could  not  locate  required  elements  from  an 
AMOS  run. 


31.  UNABLE  TO  RUN  SYSINT  FOR  THIS  CASE 

Run  of  SYSINT  VFUN  could  not  locate  elements  specified  for  detector, 
optics,  etc. 

32.  UPPER  (LOWER)  INTEGRATION  LIMIT  LOWERED  (RAISED)  TO  NNN  NM 

Run  of  SYSINT  VFUN  adjusted  integration  limits  from  400-1100 
nanometers  to  conform  to  range  of  one  or  more  of  the  specified  elements. 

In  addition  to  these  messages  from  the  AMOS  package,  many  error  messages  may 
be  generated  by  the  SPECTR  package,  particularly  if  the  user  generates  his 
own  user  function,  or  does  spectral  computations  (e.g.,  GENERATE  X = DH  * .4. 
These  messages  are  generally  self-explanatory. 

Other  error  messages  such  as"JCL  ERROR"  or  "INSUFFICIENT  CORE"  may  be  generated 
by  the  computer  operating  system,  and  will  depend  upon  the  installation. 


APPENDIX  C 


SYSINT  VFUN  DEFAULT  ELEMENTS 


The  tables  of  spectral  data  default  elements  shown  in  Tables  C-l 
through  C-3  are  called  in  SYSINT  VFUN  when  no  specific  elements  have 
been  named  in  the  command.  The  default  value  of  the  system  constant 
is  1.0.  These  spectral  data  also  illustrate  the  I/O  PRINT  XXX  command 
discussed  in  Section  6 of  this  report. 
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TABLE  C-l 

DEFAULT  SPECTRAL  ARRAY  FOR  THE  DETECTOR  (PEAK  NORMALIZED) 


VALU63  NM  IWN 


400 

0.4V30 

*70 

0.0733 

000 

0.7012 

10*0 

0.1307 

*00 

0.7031 

*70 

0.07*3 

040 

0.7*01 

10*0 

0.1217 

400 

4.71*0 

400 

0.07*3 

040 

0.7401 

1030 

0.1127 

000 

0.7271 

400 

0.0010 

070 

0.7310 

1030 

0.1044 

300 

0.73V* 

400 

0.0030 

*70 

0.7120 

1040 

0.94376-01 

010 

0.7012 

400 

0.0000 

*00 

0.4020 

10*0 

O.07436-OI 

010 

0.7*33 

700 

0.0*70 

0*0 

0.4720 

1070 

0.7*406-01 

*20 

0.770* 

700 

0.0000 

*00 

0.4000 

1070 

0.73076-01 

020 

0.700V 

710 

0.0010 

*00 

0.4276 

1000 

0.4**06-01 

000 

0.0427 

710 

0.0034 

ooo 

0.4020 

1000 

0.42406-01 

000 

0.0170 

720 

0.0034 

ooo 

0.3*04 

1000 

0.0442E-01 

040 

0.0310 

720 

0.0070 

010 

0.0300 

1000 

0.01346-01 

0*0 

0.0*14 

730 

1.000 

OlO 

0.0321 

1100 

0.44376-01 

ooo 

0.0010 

730 

0.0000 

020 

0.3030 

000 

0.0*11 

740 

o.oooo 

020 

0.40*0 

* ( 

040 

0.0701 

740 

o.oooo 

030 

0.443* 

0*0 

0.0773 

700 

0.00*0 

030 

0.4*3* 

*70 

0.00*4 

700 

0.00*7 

V40 

0.4234 

*70 

o.ovoo 

7*0 

3.0*07 

040 

0.40*0 

0*0 

0.0000 

740 

0.0000 

OOO 

0.300* 

**N» 

0.0402 

000 

4.0610 

770 

0.0700 

003 

0.3*03 

0.000* 

*40 

0.0040 

77* 

0.071* 

040 

0.3*03 

*10 

0.0000 

000 

0.V137 

7*0 

0.0*30 

0*0 

0.3272 

*10 

0.0*** 

*44 

0.0107 

70S 

o.oooo 

076 

0.300* 

WO 

0.07*0 

400 

0.02*7 

TOO 

0.V477 

070 

0.2043 

«2I 

0.0*2* 

410 

0.0307 

700 

0.0377 

0*4 

0.2703 

WO 

0.0010 

410 

0.0300 

000 

0.0277 

0*0 

0.2*02 

*30 

0.0007 

*20 

0.0400 

OOO 

0.0174 

OOO 

0.2011 

0*0 

0.40*4 

420 

0.04*1 

010 

0.0074 

ooo 

0.2340 

4*3 

0.0170 

400 

0.0472 

010 

0.00*0 

1600 

0.2211 

wo 

0.4240 

*00 

0.0010 

020 

0.0034 

1003 

0.2000 

*00 

0.4*02 

440 

0.0001 

020 

0.0703 

1010 

0.104* 

400 

0*4403 

*40 

0.0002 

030 

0.0343 

1010 

0.1023 

*05 

0.4021 

*00 

0.0412 

030 

0.0*10 

1020 

0.1*41 

470 

0.4*11 

400 

0.0447 

040 

0.0273 

1020 

0.1400 

47% 

0.4717 

4*0 

0.0*73 

043 

0.0123 

1030 

0.1000 

400 

0.4027 

440 

0.0703 

030 

0.7072 

1030 

0.1407 
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DEFAULT  SPECTRAL  ARRAY  FOR  LENS  TRANSMITTANCE 


VALUE  & FUR  4JH 


445 

0.4050 

*70 

0.4110 

095 

0.0040 

1*40 

440 

o.eiwo 

*75 

0.4000 

040 

0.0000 

10*9 

4*5 

0.0150 

400 

0.0U5O 

0*5 

o.oooo 

1050 

904 

0.4*00 

445 

0.0020 

470 

0.0000 

1009 

509 

0.4210 

4*0 

0.7V40 

475 

0.4041' 

1040 

910 

0.4*20 

4*5 

0.7*70 

044 

0.0400 

1045 

515 

0.4250 

700 

0.7*50 

409 

0.0000 

1070 

520 

0.4*40 

705 

0.7*40 

0*0 

0.0000 

1075 

929 

0.4250 

710 

0.7*50 

0*0 

0.00*0 

1000 

550 

0.4250 

719 

0.7*20 

40* 

0.4004 

1045 

555 

0.4250 

720 

0.7V 14 

*05 

0.0000 

10*0 

940 

0.4250 

725 

0.TV10 

*10 

0.4000 

10*9 

549 

0.0*50 

750 

0.7*10 

*19 

0.40*0 

1100 

550 

0.0290 

755 

0.7*20 

*20 

O.tOuO 

559 

0.4290 

740 

0.7*50 

*29 

0.0040 

540 

0.4290 

745 

0.7*40 

*50 

0.4000 

945 

w.4250 

750 

0.7*50 

*55 

0.0000 

970 

0.4250 

755 

0.7*50 

040 

0.0000 

575 

0.4250 

740 

0.74*0 

045 

0.0000 

9*0 

0.0250 

749 

0.7*40 

090 

0.4000 

400 

4.4200 

909 

0.4*50 

170 

0.7*70 

*09 

0.7**0 

400 

0.4500 

9*0 

0.4250 

775 

0.7*00 

*40 

0.7*00 

410 

0.9200 

5*5 

0.4250 

700 

0.7*00 

*45 

0.7*70 

*15 

0.5100 

*00 

0.4*50 

70S 

0.7V*0 

*74 

0.7*00 

420 

0.4100 

*05 

0.0250 

TOO 

0.79*0 

*75 

0.7*50 

420 

0.4400 

410 

0.0*50 

7*5 

0.4004 

*40 

0.7040 

400 

0.4700 

419 

0.0250 

000 

0.4000 

*05 

0.7*50 

459 

0.4*00 

420 

0.4250 

005 

0.0000 

v*o 

0.7*20 

440 

0.7*00 

425 

0.4*50 

010 

0.4*00 

**s 

0.7*10 

449 

0.7990 

450 

0.0240 

015 

0.4000 

1000 

0.7*00 

440 

0.7900 

*55 

0.4250 

020 

0.8000 

1009 

0.74*0 

459 

0.1400 

*40 

0.0220 

029 

0.4000 

1010 

0.7400 

440 

0.770* 

445 

0.0210 

050 

0.4000 

1019 

0.7470 

449 

0.7440 

•90 

0.424)0 

459 

0.4000 

1424 

0.744* 

470 

0.7440 

*59 

0.01*0 

440 

0.4004 

1025 

0.7090 

479 

0.74  50 

*40 

0.4170 

049 

0.0000 

1050 

0.7040 

400 

0.4000 

4*9 

0.4140 

090 

0.4000 

1059 

0.7050 

».T»t 
0.7410 
• •WOO 
••TWO 
0.77*0 
0.77*0 
0.7720 
0.7700 
0.7400 
0.7440 
0.7440 
0.7420 
0.7400 


TABLE  C-3 


DEFAULT  SPECTRAL  ARRAY  FOR  STANDARD  GROUND  REFLECTOR 


VALUES  FOR  025F 


AOS 

0.0*001-01 

07* 

0 

1272 

OSS 

0 

144* 

1040 

0 

1*4* 

4*0 

0. *0*01-01 

479 

0 

1277 

040 

0 

1474 

1045 

0 

1472 

4*5 

0.01*01-01 

400 

0 

1202 

049 

0 

1400 

1090 

0 

1*77 

500 

0.93001-01 

409 

• 

1200 

070 

* 

1405 

1055 

0 

14*2 

SOS 

0.04001-01 

4*0 

0 

12*3 

075 

0 

14*0 

10*0 

0 

140* 

S10 

o.vsiet-oi 

4*5 

0 

1200 

0*0 

0 

14*4 

10*5 

0 

1**3 

SIS 

0.9*101-01 

700 

0 

1304 

005 

0 

1501 

1070 

0 

1**0 

520 

0.972*1-41 

70S 

o 

130* 

0*4 

0 

1944 

1075 

0 

1704 

SSS 

0.90301-01 

710 

0 

1314 

4*5 

• 

1512 

10*0 

o 

170* 

sso 

0.V.T3  01*01 

715 

0 

1320 

*00 

o 

1917 

10*5 

o 

1714 

0.6lbbt-01 

SSS 

0.1004 

720 

* 

1325 

*05 

0 

1522 

10*0 

0 

1720 

)»$ 

0.4204t-01 

540 

0.1415 

72S 

0 

1330 

*14 

0 

1520 

10*5 

0 

1725 

360 

0.43101-01 

S4S 

0.1029 

730 

0 

1334 

*15 

0 

1533 

11O0 

0 

1730 

us 

U.*420t-01 

SSO 

0.1034 

735 

0 

1341 

*20 

o 

153* 

1105 

o 

173* 

370 

0.4S20E-OI 

SSS 

0.1047 

740 

0 

1944 

*25 

0 

1544 

1110 

0 

1741 

STS 

0.44  SOF -01 

544 

0.1057 

745 

0 

1352 

*90 

o 

154* 

1115 

0 

1744 

Mb 

0 .47 401-01 

S4S 

0.2040 

750 

0 

1357 

*35 

0 

1544 

1120 

0 

1752 

US 

0 .44441-01 

S70 

0.107* 

799 

0 

1342 

*40 

0 

1540 

1125 

0 

1757 

**o 

Cl  .4*  SOt -01 

STS 

0.100* 

740 

0 

1340 

*49 

0 

1S45 

1130 

0 

1742 

MS 

0.70404-01 

500 

0.1100 

745 

0 

1373 

*30 

0 

1370 

1135 

• 

17*4 

400 

0.71401-01 

SOS 

0.1111 

770 

0 

1ST* 

VSS 

0 

1574 

1140 

0 

1773 

AOS 

0.72704-01 

5*0 

0.1121 

775 

0 

1304 

*40 

0 

1501 

1145 

0 

1770 

AiO 

0 .73001 -01 

STS 

0.1132 

700 

0 

130* 

*49 

o 

1504 

1150 

0 

17*4 

A1S 

0.74001-01 

400 

0.114S 

705 

0 

13*4 

*70 

0 

15*2 

1155 

0 

ITS* 

420 

O.?5*bt-Ol 

403 

0.1153 

7*0 

0 

1400 

*75 

o 

15*7 

1140 

0 

17*4 

ASS 

0. 77001 -ol 

410 

0.1144 

7*5 

0 

1405 

*00 

• 

1402 

11*5 

0 

1*00 

430 

0.70001-01 

41S 

0.1175 

*00 

0 

1410 

*05 

0 

1400 

1170 

0 

1005 

ASS 

0.7*101-01 

420 

0.1105 

005 

0 

1414 

«*0 

0 

1419 

1175 

0 

1010 

440 

0 **w2wt— 41 

423 

0.11*4 

010 

0 

1421 

**s 

0 

1410 

1100 

0 

101* 

44  S 

0.01201-01 

490 

0.1207 

SIS 

0 

1424 

1000 

0.1424 

1103 

o 

1021 

4SO 

0.02301-01 

4SS 

0.1217 

020 

0 

1432 

1005 

0.142V 

11*0 

o 

1024 

ASS 

0.03401-01 

440 

0.1220 

025 

0 

1437 

ioio 

0. 

1*34 

11*3 

0 

1032 

440 

0.04401-01 

445 

0.12S* 

030 

0 

1442 

1015 

0. 

14*0 

1200 

o 

1037 

445 

O.OSSOl-Ol 

400 

0.1250 

OSS 

4 

1440 

1020 

0. 

1*45 

4T0 

0.04401-01 

OSS 

0.12S4 

040 

0 

1459 

1025 

0.1450 

ATS 

0. 04701-01 

440 

0.1241 

04S 

0 

1450 

1090 

0.1454 

400 

0.00700-01 

44S 

0.1244 

090 

0 

1444 

1035 

0.1441 
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PROGRAMMER'S  AIDS 


This  appendix  contains  guidelines  for  the  programmer  who  will  main- 
tain SCAT3.  The  call  hierarchy  of  SCAT3  subroutines  is  shown  in 
Figure  D-l,  subroutine  definitions  are  contained  in  Figure  D-2;  and 
the  JCL  to  run  program  SCAT3  is  shown  in  Figure  D-3. 


I 


D-l 


"SPECTR" 


GENERA 


SETUP 


ECRITE 


DELETE 


CFTOK 


PCHEK 


OZONE 


CKDATE 


CREARH 


ATDF.N 


WATER 


ALBRF 


TRANS 


TRANA 


Figure  D-l.  Call  Hierarchy  of  SCAT3  Subroutines 


BEMAS 


POLCOR 


BEMAS 


-2 


Figure  D-2.  Program  Notes 


SCAT3  is  a generalized  program  for  manipulating  spectral  elements,  and 
generating  optical  characteristics  of  the  atmosphere. 


SUr&OUTf NlS 

ft********** 

A.  tT V,)SP^c'r> I c PACKAGE.: 

a 'os  "v.\inline". .interface  with  spectr.  decodes  data  card 

CHECKS  6/OP  SETS  UP  CONSTANTS  6 DEFAULTS.  SETS 
Un  TEMPORARY  ELEMENTS  IN  COKE  DATA  TABLE. 


A L ^ F GENERATES  APPROPRIATE  ALBEDO  SPECTRAL  REFLECTANCE 

A SOL  CALCULATES  AEROSOL  SCATTERING  PARAMETERS  (KA«C> 

ATOff.  DETERMINES  RAYLEIC-H  SCATTERING  COEFFICIENT  <KR> 

3E*-' AS  AIR  MASS  PATH  LENGTH  CORRECTION  tSECANT  THETA) 

BLOCK  DATA  A IRMAS 

CFTGK  CONVERT  FEET  TO  KILOMETERS 

CKCATE  CHECKS  INPUT  DATA  FOR  DATE 

CREAR-.  CALCULATES  PKECIPTIBLE  MOISTURE  6 RH  WHEN 

TS5C  IS  SPECIFIED.  MAY  INCREASE  DEFAULT  TEMP 

D-<  GENERATES  SPECTRAL  ELEMENT  FOR  DAYLIGHT  I RR  AD  I ANCE 

ON  A HORIZONTAL  SURFACE 

D7  GENERATES  SPECTRAL  ELEMENT  FOR  DAYLIGHT  IRRADIANCE 

ON  A VERTICAL  SURFACE 

OZONE  ESTIMATE  OF  OZONE  BASED  ON  DATE  6 LATITUDE 
BLCCK  DATA  OZT ABL 

PCHEK  CHECKS  INPUT  PRESSURE  6 CORRECTS  FOR  ELEVATION 

polcor  computes  angular  scatter  function  <"fish"> 

SH  GENERATES  SPECTRAL  ELEMENT  FOR  SKYLI6HT 

IRRADIANCE  ON  A HORIZONTAL  SURFACE 


D-3 


Figure  0-2  (Continued) 


s<  generates  spectral  element  for  skylight  radiance  (haze> 

sve  GENERATES  SPECTRAL  ELFMENT  FOR  SKYLIGHT  IRRADIANCE 

ON  A VERTICAL  BACKLIT  SURFACE  (SHADOW) 

S'J\  GENERATES  SPECTRAL  IRRADIANCE  OP  SUN  (ABOVE  ATMOSP-iERF) 

TRANA  CALCULATES  SPECTRAL  ApSOPPTlONS  OF  WATER*  OZONE  6 OXYGEN 

TO A NS  CO  1PUTFS  SPECTRAL  TRANSMITTANCE  DUE  TO  SCATTERING 

aATER  CALCULATES  PRP C I PT ] 3LE  MOISTURE  ( MOISTURE  SCALE  HEIGHT) 

VFONCI  USER  FUNCTION  WHICH  PERFORMS  "SySINT  VFUN" 

SYSTEM  INTEGRATIONS 


~.  FRC  SPECTR: 

I \ PR  ?GR<'M  MAINLINE*. TRANSFERS  INf-UT  COMMAND  AND 
DATA  CARDS  FROM  //FT05  FILE  TO  //FTC2  FILE. 

"OPENS"  DISK  LIBRARY  FILE  //pTPB  AND  READS  DIRECTORY 
DECODES  COMMANDS  AND  DIRECTS  EXECUTION  TO  THE 
PROPER  SUBROUTINE. 

GENERA  DECODES  6 PROCESSES  GENERATE  COMMAND 

IORC  DECODES  6 PROCESSES  I/O  COMMAND 

FILE  DECODES  6 PROCESSES  FILE  COMMANDS 

NOTE:  AMOSand GENERATE  SCAT 3 AMOS  COMMANDS  ARP  DIRECTED 

TO  SUBROUTINE  A''OS 

SETUP  ENTERS  AN  ELEMENT  NAME  IN  THE  CORE  DATA  TABLE 

LDCATt  LOGICAL  FUNCTION  TO  SEARCH  FOR  AN  ELEMENT  NAME 

IN  THE  CORE  DATA  TA5LE*  OR  IN  THF  DISK  LIBRARY. 

DELETE  RF'IOVF  S AN  ELEMENT  NAME  FROM  THE  CORE  DATA  TABLE* 

BUT  NOT  FROM  THE  DISK  LIBRARY. 

ECRITE  WRITES  AN  ELEMFNT*  ARRAY  OR  CONSTANT*  ON  THE 
LINE  PRINTER. 

THERE  ARE  NUMEROUS  OTHER  SUBROUTINES  IN  SPECTR  WHICH 

support  those  already  mentioned,  consult  the  program 


LISTINGS  FOR  DETAILS. 


end 
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Figure  D-3.  JCL  to  Run  Prograa  SCAT3 


C 

//* 

//• 

/ /* 

//STE°ONF  EXEC  FuPTHC  « PAR*. FODT  = • MAP . XRE* .OPT  = 2 • 

//FJRT.SYSIN  DC  * FOR  TR  AN  SOURCE  DECKS  FOLLOW  THIS  CARD 
/ /FOR  T • SY  SUT 2 DD  UN  I T = SYSDA . SPACE = ( C YL » 20 »RLSE ) 

/* 

//STEPTWC  exec  AiMFC»PARM.ASM='LIST»LOAD' 

//ASM.SYSGO  DD  DSN=*.STEPONE.FORT.SYSLIN»DISP= (MOD. PASS* DELETE) 
//ASV.SYSU  DD  * ASSEMBLER  SOURCE  FOR  "CORE"  FOLLOWS  THIS  CARD 

/* 

//STEP3  EXEC  Ai"FC.PAR«.ASM='LIST.LOAD' 

//ASV.SYSGO  DC  DoN=*.STEPONE.RORT.SYSLIN,DISP= ( MOD » PASS ♦ DELETE > 
//AS”.SYSI\  DC  * ASSEMBLER  SOURCE  FOR  "DATE"  FOLLOWS  THIS  CARD 
/* 

//STEPFOUR  EXFC  FORTHLG.PaR*.LKED= (XREF.LIST.LET) 

//LKE^.SYSLIU  DD  DSN  = * . STF.PONE  . FORT  . SYSL  I N »D  I SP=  ( OLD  .DELETE  ♦ DELETE ) 
//  DD  * LINK-EDIT  INPUT  FOLLOWS  THIS  CARD 

l\TRY  main 

/* 

//G0.FT02F00I  DD  UNI T = SY SDA » SPACE = (CYL .( 10 » 1 >> . 

//  DC?  = (RECFM=F-?.LRECL  = 8C.5LKSIZE  = 3600) 

//GO.FTORFOOl  DD  DSN= SPECTR AL . DAT A .UN I T =DI SK » VOL= SER= *BF225 . D I SP=OLP 
//GO.FllCFOOl  DD  UN  I T = SYSDA  ♦ SPACE  =(  732  «(  450  * 1 > > 

//GO. SYS  I N DD  * INPUT  DATA  CARDS  FOLLOW  THIS  CARD 

/* 

ALL  JCL  CONVENTIONS  ARE  FOR  THIS  INSTALLATION'S  IBM  370.  RUNNING 
UNDER  OS.  USING  THE  LASP  OPTION. 

PLEASE  NOTE  ThE  FOLLOWING..... 

/*S£TUU • • • 

TrIS  CARD  PAY  3E  REQUIRED  FOP  OPERATOR  MOUNTING  OF  THE  VOLUME 
IN  THE  'ID'  PARAMTER*  WHICH  IS  A DISK  PACK  ON  WHICH 
RESIDES  THE  SPECTRAL  DATA  LIBRARY  (SEE  FT09F001  OF  THE  'GO* 
STEP) 

//STEPONt ... 

SEVFPAL  OF  THE  FORTRAN  ROUTINES  TAKE  ADVANTAGE  OF  THE  POWER 
OF  THE  I EM  'H'  LEVEL  COMPILER*  AND  WILL  NOT  COMPILE  UNDER 
A LOWER  LEVEL « WITHOUT  MODI  FACTION.  THE  COMPILER  OPTIONS 
STATED  IN  THE  'PaRM'  FIELD  ARE  HIGHLY  RECOMMENDED. 

TnlS  JCL  IS  REQUIRED  FOR  OUR  CATALOGUED  PROCEDURE  INVOKING  THE 
• H • COMDILEK. 

//steptwp... 

THIS  IS  OUR  CATALOGUED  PROCEDURE  TO  INVOKE  The  'IEUASM' 

LANGUAGE  COMPILER.  NOTE  THE  ADDITION  OF  THE  PROCESSOR'S  OUTPUT 
TO  THAT  ALRtADY  PRODUCED  BY  THE  FORTRAN  COMPILER. 

//STEP3. • . . 

SAME  AS  EXPLAINED  FOR  STEPTWO  (ABOVE) 

//STEPFOUR.... 

THIS  IS  OUR  CATALOGUED  PROCEDURE  TO  INVOKE.  THE  PROGRAM  • EWL  • 
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Figure  D-3  (Continued) 


(IP-V  LINKAGE  EDITOR).  DEPENDING  UPON  THE  PARTICULAR  SYSTEM 
being  USED»  AN  OVERLAY  STRUCTURE  may  be  advisable  TO  CUT 

DOWN  THE  ME  "OR Y REQUIREMENTS.  HOWEVER , IF  THE  PROGRAM  CAN  EE 
EXECUTED  AS  A SINGLE  SEGMENT « RUN  TIME  WILL  BE  OPTIMIZED. 
//GO.CTOZFOOI... 

TEMPORARY  DATA  SET,  NO  NEED  TO  ALLOCATE  PERMANENT  SPACE  FOR 
IT.  THIS  SET-UP  SPECIFIES  THAT  THE  DATA  SET  IS  NEW,  AND  WILL 
BE  SCRATCHED  AFTER  THIS  RJN. 

//GQ.FTC9F001.... 

THIS  IS  THE  SPECTRAL  DATA  LIBRARY,  WHICH  WF  KAVF  SET  UP  ON 
A 3330  DISK  PAC*,  VOLUME  SERIAL  NUMBER  ABF 225 

THIS  PARTICULAR  SETUP  ASSUMES  THAT  THE  DATA  SET  HAS  BEEN  PRE-FOR *ATTEC 
O'.  The  DEVICE,  WITH  A SPACE  PARAMETER  OF  (7a0,Z920). 

THE  PRE-FOR "ATT  I NG  COULD  HAVE  BEEN  DONE  BY  A PREVIOUS  RUN, 

WITH  DlSPsNEW,  OR  WITH  • I EHPROGM * UTILITY  PRCGR AM • 

//G0.FT1CF001.... 

This  is  ANOTHER  SCRATCH  DATA  SET  USED  EY  THE  PROGRAM. 


SINCE  THE  PROGRAM  IS  QUITE  LARGE , IT  IS  HIGHLY  ADVISED  THAT 
A LOAD  MODULE  PE  CREATED,  AND  RFTAINED  IN  A PRIVATE  LIBRARY, 

AS  SOON  AS  POSSIBLE. 

NOTE  THE  COUNTS  IN  THE  FORTRAN  MAINLINE,  REGARDING  INITIAL 
USE  OF  the  spectral  element  LIBRARY  ( preceeding  THE  CALL  TO 
•FDTI.NT*  ) . 

THIS  PROGRAM  SHOULD  RUN  ON  ANY  INSTALLATION  SUPPORTING 
IBM  OS/ 360-370 , ASSUMING  ENOUGH  MAIN  STORAGE,  HAVING  SUITABLE 
RANDOM- AC CESS  STORAGE  FOR  THE  DATA  LIBRARY,  WITHOUT  ANY  SOURCE 
MODIFICATION.  THE  JCL  CAN  EASILY  BE  MODIFIED  TO  ACCOMODATE  ANY 
VERSION  OF  ASP,  OP  a SYSTEM  RUNNING  WITHOUT  THE  BENEFIT  OF  ASP. 
C 


— 
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APPENDIX  E 

USER  WRITTEN  FUNCTIONS 


The  user  may  wish  to  write  special  functions  to  perform  repetitive  calcu- 
lations on  elements  generated  by  SCAT3,  or  provide  outputs  tailored  to  a 
particular  job.  These  may  be  written  as  FORTRAN  SUBROUTINE  UFUNCT,  VFUNCT, 
or  WFUNCT.  These  SUBROUTINES  may  be  invoked  by  a command  card  of  the  form: 

GENERATE  XXX  UFUN  YYY  ZZZ  AAA=N 

where  XXX,  YYY,  ZZZ,  and  AAA  are  arguments  that  might  be  passed  to  SUBROUTINE 
UFUNCT.  Similarly,  GENERATE  XXX  VFUN  will  result  in  a CALL  to  SUBROUTINE 
VFUNCT,  and  GENERATE  XXX  WFUN  will  invoke  SUBROUTINE  UFUNCT. 

SCAT3  makes  use  of  the  above  capability  in  the  GENERATE  SYSINT  VFUN  command. 

It  is  recommended  that  before  one  begins  writing  a user  function,  a FORTRAN 
programmer  review  the  program  listings,  particularly  SUBROUTINE  VFUNCT  and 
SUBROUTINE  AMOS.  It  should  be  noted  that  the  SPECTR  portion  of  the  SCAT3 
code  contains  dummy  versions  of  SUBROUTINE  AMOS,  UFUNCT,  VFUNCT,  and  WFUNCT. 
For  this  reason,  when  one  performs  a LINKAGE  EDIT  to  incorporate  a new  user 
function,  the  new  SUBROUTINES  must  be  loaded  first.  Numerous  utility  SUB- 
ROUTINES within  the  SPECTR  coding  may  be  used  for  arithmetic  and  bookkeeping 
functions.  At  such  times,  the  listing  for  details  and  calling  sequences 
should  be  consulted. 

As  a guide  in  writing  a user  function,  a simple  card-punching  aid  is  described 
below.  A listing  of  the  FORTRAN  statements  can  be  found  in  Figure  E-l. 

The  user  function  can  be  used  to  punch  any  spectral  array  element  previously 
calculated  or  that  is  in  the  library,  onto  cards  in  the  format  used  to  enter 
a new  element  via  an  I/O  READIN  command.  The  element  values  will  also  be 
printed.  The  function  will  be  invoked  by  the  command: 

GENERATE  XXX  WFUN  FORM=N. 

where  XXX  is  the  name  of  the  array  element  (e.  g. , the  AMOS  output  array  code  DH) 
and  N is  a format  code.  In  addition  to  illustrating  the  writing  of  a user 


E-l 


function,  the  example  may  be  of  practical  use. 


The  following  steps  are  performed  in  this  user  function: 

(a)  Check  the  number  of  arguments. 

(b)  Use  SPECTR  function  LOCATE  to  determine  the  position  of 
the  values  for  element  XXX  in  COMMON/ L I BRY/ARRAY(.)  if  the 
element  name  is  not  in  the  core  data  table,  the  disk 
library  is  searched  and  the  element  loaded  into  ARRAY. 

(c)  If  the  element  cannot  be  located  or  if  the  mode  is  wrong; 
e.g.,  it  has  been  defined  as  a constant  instead  of  an 
array,  an  error  message  is  generated. 

(d)  The  LOW  and  HIGH  positions  in  the  element  are  determined, 
and  the  date  is  obtained  from  the  system. 

(e)  Check  the  parameter  "FORM  *"  for  1.,  2.,  or  3.  A code  of 

1 will  use  FORMAT  (10G7.2)  and  a code  of  2 will  use  FORMAT 
(10E7.2).  Code  3 will  use  FORMAT (10F7. 4) , which  may  not  be 
sufficient  for  some  element  values. 

(f)  Punch  the  array  values. 

(g)  Print  the  array  values. 

The  SUBROUTINE  may  be  compiled  and  linkage-edited  to  SCAT 3 . The  new  sub- 
routine must  be  loaded  before  SC AT 3,  which  contains  a dummy  routine  by  that 
name.  The  JCL  must  be  modified  to  include  a card  punch  as  device  FT07F001. 


15*  000? 


r$*1>0T3 

IS*  000* 

- is*  'oeo5 

lS*_®80f> 
TS*  0007 

— l wo  oar 

is*  OOOO 


~ 7SN~0010 

— rwTtn?' 

IS*  0013 
— ISNT1D1* 

— IT*  "6  01  S' 


Figure  E-l 

FORTRAN  Listing  of  Sample  User  Function 

SUBROUTINE  WFjnCT  (RESUlT. PARS, NPAP) 


USE*  FUNCTION  TO  PUN C-  CARDS  FRO*  ARRAY  ELEMENT 


) 1067.3 

10ET.3 

1 OFT.i  IUSE  wITm  CASE.. EASILY  OVERFLOWED  1 ) 


LOGIC*:.  LOC»Tl 

RtAL*f>~l>ESULTiPAAS  (*AaB)  .TOURS 

CT*f  NSTCsToTT  CT2) 

CQhmqn/l IBRY/AR R*v (163. ISO) » CONST  1600) 
FOuNalEnce  (array  (i,n  .range  (1 . 1 ) ) 
CO*MO*/I*PUT/STRj*r.  (801  .lSTR.INFI 


• OF  arguments  ok  ? 

rrcN&AR  .Tf  .ITStTTo-! 


on  F6RpaT(«D*.1Da.**»* — T 00  * any  ACUMEN 
WRITF (S.lol) 

— - n^roRk,- 

C 6ET  FLEMENT  with  name  'MRFSUlT)" 

1 Ir  (LOCATE  (T»E5ulT  . JPOSTmOBe  ) ) GO  to  ? 


SUMENTS  ) •/•('«) 


I5N  0016  WRITE tS.10?)«ESULT 

i s*  goto urtuRR 

C MODE  OK  r 


LOCATE  • * AP ) 


IS*  0023 
— 1SN  062a’ 

fSN~b®25~ 
1SN  OQgfe 

IS*  0027 
— ISN  8026~ 

ISN  002B 

1 5*  0030'- 
15N  0031 

ISN  0032 
— lSN~0033 

“I5H“®03T- 
1SN  003» 
ISN  TOOT7 


< »0*. 10*. ••••  * « A6 « ' NOT  AN  ARRAY  ElEMENT»/»0* > 

WRITE (6.103) 

urTjRk, : 

_C SET  LOW  L hIGh  INDEX  FOR  ELEMENT  ILAST  2 ARGUMENTS) 

5 l6w*banGE (162.JP0S) 

IH1GH«RANC.E  ( 163.  JPOS) 

t "cOnVeAT” to  wavElEnOIm:  U)  ■ 300  NANOMETERS 

LAML«2»S  » 5*L0W 

sh 

C GET  DATE  FROM  SYSTEM  VIA  SPECTR  UTILITY 

nirtofiomionn 

C CHECK  FORMAT  OPTION 

irOBM.l 

CALL  SCALNM(STRING.S0.'  FOR*. • ,6,F ,i« ) 

C NDTHOE  ACRTD  FOUND 

1F0PM» IFI R <F* . 1 ) 

Tf I VlFdffMTbt • 1 ) • AND. < 1FORM.LC • 31 ) 60  TO  A 

c 

WRITE (S. 16411 

104  FORMAT(*0».10A,**w*  INVALID  FORMAT  COOE'.FO.l.f  <1.-3,  OKIl/IO*) 
BFTUBIf 
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Figure  E-l  (Continued) 


rr— i 


tSN  003B 
ISN  6639 


TSN  00*0 
"YSN  6041  ' 
ISN  004? 
ISN  0043 
_ISN  004* 

TSN  004S 


ISN  00*0 
ISN  06*7 
ISN  00*0 
'ISN  0049 
ISN  0050 

Twwri — 

ISN  005? 
ISN  0053 


1 *Wm"«E4&E* 

4 «BI TE ( T« 1 05 ) RESULT. L»NL*L*NN.IO»TE 

Yo5  ?0P*»A7<»I/0  PET6FR*/ 

1 40 • • V4l*  *0X  * 1 4 • 0 X * I4*40X«?44 ) 


SET 

FID«PESULT 

“Rf 4i*0  • t Ihigh-lO«i /l6  * 1 

I 1*L0p 

'DO  6 N.r.NCiSD 
I?«  U »9 

•*i&rroFBm 


60  T0(S.b,7) t Ir09M 

T 

100  «'OP*'4T(1067.?t?*»4*tI*) 

5 *e i tftt, io-sni»wT(TTj*Rssr,T«i i .ur.no. n 

60  TO  8 

'107  r On k 47 TT6 £7 .?•  ?T;'4*TT*1 

0 «OITE (7,1071 (40P4Y ( I . JPOS1 .1*11.1?) .PIO*N 

ocm’n 


100  »’OP-»T(10P7.*.?»,44,I4) 

7 ppite (7,i68nApniir(r.jP05rryiniti?T.TTCi;v 


ISN  005* 


c pn i nt 

ISN  0055  4PITE(0.1O91PE5UL7 

isn  ooss  roo  eePMATTn'/*(rr,ior,>pjKCHEo  T7i*7‘irn 

C •■PESulT"  » I lL  9E  OPINTEO  UPON  PETUPN 

ISN  0057  netoPN 

I S N 00  50  EJ*0 


